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ABSTRACT 

Numerical simulations of cosmological structure formation show that the Universe’s most massive 
clusters, and the galaxies living in those clusters, assemble rapidly at early times (2.5 < z < 4). 

While more than twenty proto-clusters have been observed at z 2 based on associations of 5-40 
galaxies around rare sources, the observational evidence for rapid cluster formation is weak. Here we 
report observations of an asymmetric, filamentary structure at z = 2.47 containing seven starbursting, 
submillimeter-luminous galaxies and five additional AGN within a comoving volume of 15000 Mpc 3 . 

As the expected lifetime of both the luminous AGN and starburst phase of a galaxy is ~100Myr, 
we conclude that these sources were likely triggered in rapid succession by environmental factors or, 
alternatively, the duration of these cosmologically rare phenomena is much longer than prior direct 
measurements suggest. The stellar mass already built up in the structure is ^10 12 M 0 and we estimate 
that the cluster mass will exceed that of the Coma supercluster at z ~ 0. The filamentary structure 
is in line with hierarchical growth simulations which predict that the peak of cluster activity occurs 
rapidly at z > 2. 

Subject headings: galaxies: clusters: general — galaxies: starburst — galaxies: quasars: general — 
cosmology: large-scale structure of universe 


1. INTRODUCTION 

An outstanding question on the study of massive 
galaxy clusters in the Universe is how and when the mem¬ 
ber galaxies built most of their mass: in a cascade before 
the cluster coalesces, gradually as the cluster accretes 
mass, or predominantly after the formation of the clus¬ 
ter. While galaxies in nearby coalesced clusters have sup¬ 
pressed star-formati on rates due to “cluster quenching,” 
(Cooper et al.|2008) it is un clear whether or n ot a rever¬ 
sal occurs at nigh-redshift (Elbaz et al. 2007), whereby 
galaxies in more massive structures have enhanced star- 
formation rates, in line wi th expectation from hierarchi¬ 
cal growth formation (Moster et al. 2013). Dense envi¬ 
ronments that are undergoing a rapid formation in the 
form of proto-clusters are difficult to detect because the 
intracluster medium at z > 2.5 has not yet been heated 
sufficiently to emit in the X-rays or absorb cosmic mi¬ 
crowave background photons via the Sunyaev-Zel’dovich 
effect. The existing discoveries of high-redshift dense 
structures do not provide adequate observational evi¬ 
dence to interpret how and when the galaxies in those 
regions formed. 

Here we present data on a distant proto-cluster at 
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z = 2.47 found serendipitously during a redshift survey 
of dusty star-forming galaxies (DSFG s) in t h e Scuba-2 - 
imaged portion of the COSMOS field (Casey et al. 2013). 
We discuss data and observations in § present relevant 
results and calculations in § [3j and discuss the implica¬ 
tions on the formation of early proto-clusters in § [4] 

2. DATA & OBSERVATIONS 

This structure, which we call PCL1002, was found 
serendipitously in Keck MOSFIRE (21-Dec-2012, 31- 
Dec-2013 and 19-Jan-2014) spectroscopic follow-up of 
ScuBA-2-selected DSFGs in the COSMOS field. The 
DSFGs’ FIR-photometry is given in Table |Tj Ob¬ 
serving conditions for MOSFIRE nights were favor¬ 
able, with clear skies and 0.5-0.7" seeing. Six DS¬ 
FGs were spectroscopically confirmed with Ha redshifts 
within ±0.007 of z = 2.472. The DSFGs’ near-infrared 
counterparts are very secure for 4/6 D SFGs, driven 
by precise 450/im positions (all given in Casey et al. 
[25l3t except DSFG J100026.73+022411.3 wlth - ^ = 
14.6 ± 4.1mJy). Two others have more ambiguity due 
to 850/im-selection. DSFGJ100018.17±022250.4 resem¬ 
bles a major-merger spanning 2" with multiple knots. 
DSFGJ100027.14±023140.8 has both 24/rm/radio emis¬ 
sion overlapping with IRAC emission towards the iden¬ 
tified counterpart, given by source name. Two other 
MOSFIRE targets were confirmed in this same redshift 
interval. The one-dimensional Ha spectra are shown in 
Figure [lj 

One additional i/ersc/iel-SPIRE-detected galaxy, 
COLDz J100018.21+ 023456.7, sits at z = 2.4790 

(Lentati et alj20151, confirmed via detection of CO(l-O) 
in a 6.5 arcmin 1 2 blank-field CO search program north of 
PCL1002 (D. Riechers etal., in prep). 

Supplementary data are pooled from a repository of 
legacy ancillary data in COSMOS. An additional 34 
spectroscopically-confirmed sources in the ^COSMO S 
survey are within 2.463 < z < 2.487 (Lilly et al. 2009), 
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TABLE 1 

Deboosted FIR-Photometric data for PCL1002’s dusty starbursts 


Name 

S 2 4 

Sioo 

Sl60 

S250 S350 S450 S500 

■5850 

Si.4GHz 

Ar 


[+Jy] 

[mJy] 

[mJy] 

[mJy] [mJy] [mJy] [mJy] 

[mJy] 

[/Ay] 

[Lq] 

DSFG J100036.03+022151.1... 

194+17 

- 

- 

11.3+2.2 15.6+2.7 11.3+4.9 14.5+3.1 

4.6+1.1 

- 

(3.12+^)xl0 12 

DSFG J100018.17+022250.4... 

128+16 

- 

- 

11.3+2.2 - 3.9+4.1 

3.3+1.0 

- 

(3.84l?i®)xl0 12 

DSFG J100016.57+022638.4... 

890+17 

6.7+1.9 

19.2+3.6 24.5+2.2 21.5+2.7 17.3+4.7 11.6+3.0 

3.7+1.0 

5716+73 (7.52+?'H[)xl0 12 

DSFG J100056.83+022013.3... 

90+27 

- 

- 

9.1+12.9 18.3+6.0 

10.9+1.1 

- 

(2.06l°g)xl0 12 

DSFG J100026.73+022411.3... 

84+14 

- 

- 

11.4+2.2 14.8+2.7 10.3+5.0 17.8+3.0 

0.4+0.8 

- 

(2.99l[fg)xl0 12 

DSFG J100018.21+023456.7... 

153+11 

- 

- 

16.6+2.0 19.5+3.4 12.0+8.8 11.2+4.2 2.57+1.74 

46+10 

(4.23l 3 '^)xl0 12 

DSFG J100027.14+023140.8... 

421+152 13.6+1.8 24.4+3.6 36.9+2.2 30.8+2.8 3.9+5.0 17.0+3.4 

5.8+1.4 

67+12 

(1-I5t^)xl0 13 


Observed Wavelength (urn) 

2.0 2.1 2.2 2.3 2.4 


DSFGJ100027.13+023138.5 z=2.483 




F DSFGJ100026J73+022411.3 z=2.477 


- J MyjjAvVvyl|\/yv^^ 


L DSFGJ100016-57+022638.4 z=2.473 
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Fig. 1. — Left: extracted one-dimensional K-band MOSFIRE spectra for proto -cluster members. Right: example spectra of PCL1002 
member galaxies in the rest-frame ultraviolet from VLT VIMOS (Lilly et al. 2009). 


identified via Lyct emission, Fell, Sin and Cil absorp¬ 
tion, also shown in Figure |T] We also make use of the 
30+ ph otometric bands of imaging data available in the 
field (|Ilbert et al.||2013h . We also dr aw on the COSMO S 
Chandra X-ray (J.5-lU keV catalog (Civano et al. 2012), 
radio 1.4 GHz catalog (Schinnerer et al. 2007), and Her- 
schel PEP/PACS and HerMES/SPIRE 100-500/im cat¬ 
alogs (Lee et al. 2013:. 

We calculate the significance of this over -density using 
methods used for other z Xi 2 structures (Steidel et al. 


1998 Chapman et al.|[2009 ) by computing the likelihood 


of observing 7 DSFGs within a A z = 0.02 interval with 
the Erlang distribution function (Eadie et al. 1971]), 

p(Az\NX) = A(AA z) N ~ 2 exp(—AA z)/(N - 2)! (1) 


where p(Az\NX) is the probability that N galaxies span¬ 
ning a redshift range A z is drawn by chance. The expec¬ 
tation for the density of galaxies per unit redshift interval 
is described by A and the calculation of the probability 
assumes no clustering as the null hypothesis. Assum¬ 
ing a volume-density of 7xlO -5 Mpc -3 for > 10 12 ' 3 Lq 
DSFG s at z ~ 2.5 (taken from the best to-date lu minos- 
ity functions of DSFGs; [Casey, Narayanan, fc C ooray| 


20141, we infer that the number of DSFGs expected 
in a redshift slice of Az = 0.02 within a 150arcmin 2 
box is A = 0.64 (in a volume ~10000Mpc 3 ). This im¬ 
plies the probability of observing seven DSFGs in this 
interval is 0.002%. This corresponds to a DSFG over¬ 
density of <5dsfg = (7 — 0.64)/0.64 = 10. For com- 
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ume (< 5 dsfg ~ 100). With only 0.6 DSFGs expected in 
the given volume, this over-density only has a 0.0001% 
chance of occurring by chance, implying that PCL1002 
is unlikely to be an artifact of incompleteness or survey 
bias. 

A similar peak is seen in the redshift distribution of 
zCOSMOS Lyman-Break Galaxies (LBGs) with a maxi¬ 
mum 5 g = 3.3. Although less pronounced than the DSFG 
over-density, the over-density of Lyman Break Galaxies 
are also spatially clustered on the sky in a fila ment ary 


structure . We use a friends-of-friends (FOF; Huchra 


& Geller 1982) algorithm to formally determine which 


LBGs are in fact members of the DSFG over-density by 
selecting sources which are within 2 Mpc proper (pro¬ 
jected on the sky and line-of-sight distance) of the DS¬ 
FGs or their immediate neighbors, or within 3 Mpc of the 
DSFGs. This is similar to low -z cluster member identifi¬ 
cation techniques but with longer ‘linking length’ accom- 


modating high_-z non-virialized structures (Chiang et al. 
2013). Table|2] lists PCL1002 members. 


FCL1002 sweeps out an effective area of 200arcmin 2 
(distributed over an area extended 25'x25') and 28 Mpc 
comoving along the line of sight, and a volume 
^15000 Mpc 3 comoving (or 400 Mpc 3 proper). The range 
of redshifts A z = 0.0239 translates to a total end-to-end 
line-of-sight velocity range of Av = 2080kms _1 . Fig¬ 
ure [2] map s the structure. Note that another over-de nse 
structure (Chiang et al. 2014, 2015 Diener et al.|2015 ) at 
z = 2.44—2.45 sits nearby (at distances of ~50-100 Mpc). 

This structure also hosts an over-density of luminous 
active galactic nuclei (AGN), an important signature 
of accreting supermassive black holes. Four galaxies 
(9.5%) are luminou s X-ray sources in th is structure (at 
Lx >10 43 7 ergs -1 ; Civano et al. 2012), a factor of 21 
higher than the expected volume "density o f AGN of sim¬ 
ilarly high luminosities (Silverman et al. 2008). Even 


among the seven DSFGs, four exhibit unequivocal AGN 
characteristics in either the X-ray, radio, optical or ul¬ 
traviolet (57%): a fraction nearly twice as high as ex- 
pecte d from previous w ork on AGN in non-clustered DS¬ 
FGs (Alexander et al.|2005|. The depth of the Chandra- 
COSMOS observations is only sufficient to detect the 
most X-ray luminous AGN, so we stack the undetected 
population to search for possible low-luminosity AGN, 
but found no X-ray emission. 

One of the seven DSFGs and one member LBG 
host 


seven 
radio-loud AGN 


(Schinnerer 
has 


et 
radio 


al. 


2007). 


luminosity 


nearly analogous to 


DSFG J100016.57+022638.4 
of Li 78 MHz = 1-9 x 10 27 WHz -1 
a local Fanaroff-Riley class II AGN (|Fanaroff & Riley 
1974|. It has the potential to become the brightest 
cluster galaxy (BCG) when the structure matures to 
z ~ 0, consistent with other high-z over-densities th at 
host single radio-loud quasars ( Venemans et al.|[2007|) . 

Beyond X-ray and radio signatures of AGN " we find 
that three DSFGs and one LBG either have very broad 
Ha lines or high [Nil]/Ha ratios consistent with AGN. 
In summary, this structure hosts 6 luminous AGN, 9 


galaxies with strong signatures of AGN, 7 submillimeter- 
luminous DSFGs, for a total of 12 exceptionally rare 
galaxies. 

3. RESULTS 

3.1. Control Galaxies outside of PCL1002. 

To draw comparisons between PCL1002’s mem¬ 
ber galaxies and field galaxies, we define a set of 
spectroscopically-confirmed galaxies that lack close phys¬ 
ical associations. The redshift range of the control sam¬ 
ple is restricted to ±200 Myr of the observed structure 
(2.30 < z < 2.66) to prevent confusion between redshift 
evolution and environmentally-driven differences. Of the 
1072 zCOSMOS sources which satisfy 2.30 < z < 2.66, 
we remove sources within 20 Mpc of PCL1002 (a fac¬ 
tor of two in comoving distance beyond the boundary of 
the structure itself) and any sources with more than two 
neighbors within 5 Mpc. The resulting sample consists of 
401 galaxies; the vast majority of these are selected via 
optical-color and appear to only differ with the PCL1002 
LBGs by environment. 

3.2. Stellar Masses & Star-Formation Rates. 

We use th e MAGPHYS spectra l energy distribution 
(SED) code ( da Cunha et al.|2008 ) to constrain the UV- 
through-far-infrared SED empirically using energy bal¬ 
ance techniques. We use a stellar synthesis template 
library as input, and attenuation is determined from a 
mix of hot/cool/PAH dust grains. Our interest in using 
MAGPHYS is threefold: to measure star-formation rates 
for galaxies not directly detected in the far-infrared, to 
estimate stellar masses for all galaxies (comparing stellar 
mass estimates of DSFGs with other techniques), and to 
compare physical properties of PCL1002 members with 
field galaxies. One galaxy, QSO J100021.96±022356.7, 
lacks a SFR estimate due to lack of far-infrared detec¬ 
tion or disambiguation of quasar-dominated optical SED 
from stellar emission. 

Star formation rates for DSFGs are measured by fit¬ 
ting simple modified blac kbodies (p lus mid-infrared pow- 
erlaws, /3 = 1.8, a = 2.0; Casey|2012| ) to all available far- 
infrared/submillimeter photometry ( Spitzer , Herschel , 
Scuba-2). We integrate under the SED from 8- 
1000 /im to get the total infrared luminosity, Lir and 
convert to a star-format ion rate using SFR/MQyr _1 = 
9.5xlO -11 Lir/Lq ( Kennicuttl 1998) adjust ed for a 
Chabrier initial mass function (|(Jhabrier||2003|) . 

LBGs in PCL1002 are more massive in stars than con¬ 
trol LBGs by a factor of 1.5±0.3. Given the uncer¬ 
tain star-formation histories of DSFGs, we check that 
the MAGPHYS-generated stellar mass estimates are not 
systematically biased by computing a rest-frame H-band 
magnitude for each DSFG, removing any mid-infrared 
dusty power-law (on average, ~50% contribution), and 
converting to stellar mass using a range of appro priate 
mass-to-light ratios (Lr/M* = 7.9+2® L 0 /mag; Hain- 
Two DSFGs lack stellar mass esti- 


line et al. 

2011 

mates; DS 

FGJ1 


mid-infrared powerlaw (with no constraint on the under¬ 
lying stellar emission) and DSFG J100027.14±023140.8 is 
blended in the near-infrared with several nearby sources. 
For the remainder, we find consistency between the 
MAGPHYS-derived stellar masses and the H-band de¬ 
rived stellar masses. 
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Fig. 2.— Three spatial projections of the z = 2.47 PCL1002 proto-cluster: on the sky (panel A), and two projections with redshift 
on the x-axis and either deck (panel B) or R.A. (panel C) on the y-axis. Thick black points are confirmed PCL1002 members, and 
small gray points are spectroscopically-confirmed galaxies in 2.30 < z < 2.66. Thin gray lines connect the two nearest-neighbors. DSFGs 
have blue boxes, radio-detected galaxies have red diamonds, and X-ray-detected gal axies have yellow triangles. Peculiar velocities are not 
constrained. Four candidate group positions are marked with large magenta crosses jDiener et al.[2013|k Panel D shows the distribution of 
spectroscopic redshifts for galaxies in within 20 / of PCL1002’s center. The over-density oi dusty slarburst, submillimeter galaxies (DSFGs; 
blue) is a factor of 10 above what is expected (dashed blue line) given the known volume density of DSFGs in the field at these redshifts. 
The observed number of LBGs (dashed black line) is a factor of 3.3 above expectation at z = 2.47. 


3.3. Submillimeter Stacking. 

To investigate low-level submillimeter emission in the 
proto-cluster’s LBGs, we stack our 450pm and 850pm 
submillinreter maps at the positions of known galaxies, 
both in PCL1002 and in the control sample. Stack¬ 
ing analysis in the submillinreter requires th e removal o f 
brigh t, significantly-detected point sources (Webb et al. 
2003) so that real DSFGs do not bias the measure¬ 


ment. Thirty-one of 34 LBG members fall within the 
sensitive area of the submillimeter maps (71 of 401 con¬ 
trol sample LBGs). Flux densities are measured by 
inverse-varia nce weighting (Viero et al. 2013 Coppin 
et al. 2015). At 850pm, we measure a flux density 
of 5^0 = 0.25 ± 0.16 mJy for PCL1002 LBGs, and 
0.11±0.13mJy for control LBGs. In other words LBGs 
in PCL1002 are brighter at 850/inr by a factor of 2.3±3.0 
(although consistent with equal flux density, the likeli¬ 
hood of greater 850/im in PCL1002 is 76%). At 450pm, 
we measure a flux density of S 450 = 0.41 ± 0.85 nrJy for 
LBGs in PCL1002, and 1.66±0.69nrJy for control LBGs. 

Although low signal-to-noise due to the small num¬ 
ber of coadded sources, these measurements together 
are suggestive that the mass of cold dust and inter¬ 


stellar medium (ISM) is potentially higher for galaxies 
in the dense structure, despite their comparable star- 
formation rates. Since galaxies’ ISM masses scale di¬ 
rectly to the_ir gas masses (wit h roughly constant dust-to¬ 
gas ratio; Scoville et alj2014), we deduce that the molec¬ 
ular gas reservoirs of the structure’s galaxies are probably 
deeper, thus their potential for heightened star-formation 
relative to similar galaxies living outside of it. Follow-up 
molecular gas measurements are needed to confirm this 
intriguing lead. 


3.4. Rest-frame Optical Morphologies. 

With Hubble Space Telescope (HST) H-band imaging 
available for 21/42 PCL1002 members (and 25/401 con¬ 
trol galaxies), we investigate morphology and interaction 



able) and an interaction class (merger, interacting pair, 
or non-interacting). Full details of both morphology and 
interaction class for member galaxies are given Table [2] 
and Figure [3j Before visual classification was carried 
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TABLE 2 

Physical Characteristics of PCL1002 Members 


Name 

z 

Type 

Luv 

SFR 

M* 

Morph 

AGN 




(L©) 

(Meyr 1 ) 

(M 0 ) 

Class 

INDICATOR 

LBG J100013.62+022604.9 

2.463 

UV 

(4.26++) xlO 10 

6.67++ 

(1.14++)xl0 19 

Disk/Int 

- 

LBG J100033.33+022159.9 

2.463 

UV 

(4.36/g'gg) X 10 10 

49 5 +0 - 6 

(4.65+??) xlO 10 

Sph 

- 

LBG J100036.90+022213.8 

2.463 

Lya/UV 

(2.24++) x 10 10 

0 1 7+O.2O 
o,±, -1.02 

(5.74+??) xlO 8 

Sph 

- 

LBG J100038.35+022216.4 

2.463 

UV 

(5.86/g'gg)Xl0 10 

25.9/^g 1 

(7.92+??) xlO 9 

Sph/lnt 

- 

DSFG J100036.03+022151.1 

2.465 

Ha 

<3.Ox 10 10 

296/gf 

(1.0ll?'??)xl0 11 

Disk/Int 

X-ray 

LBG J100018.18+022837.7 

2.466 

UV 

(4.28/g gg)Xl0 10 

12 7 +3 - 4 

1 - 2.7 

(2.37++) xlO 9 

Sph/lnt 

- 

LBG J100024.36+022236.3 

2.466 

UV 

(4.88/°™) xio 10 

25-5±|;i 

(8.79+?|) xlO 9 

Sph/lnt 

- 

LBG J100050.73+021922.4 

2.466 

UV 

(6.61++)xl0 10 

98.2++ 

(1.44+??) xlO 10 

- 

- 

LBG J100031.14+023103.3 

2.467 

Ha/Lya 

(2.6l/g'gg)Xl0 10 

0 91 +1-96 

(3.22++) xlO 10 

Sph 

- 

LBG J100051.16+022305.1 

2.467 

Lya 

(2.86+' ??)xl0 10 

17.8+'? 

(2.01+ +) xlO 9 

- 

- 

LBG J100058.80+022032.4 

2.467 

UV 

(4.90+?!)xl0 10 

63 1+ 12 
DO - -8.2 

(5.94++) x 10 9 

- 

- 

LBG J100111.03+022043.4 

2.467 

UV 

(5.69+??) XlO 10 

45 0 +21 '° 

^o.u_ 17 Q 

(1.65++)xl0 10 

- 

- 

LBG J100100.91+021927.3 

2.469 

UV 

(5.07+??) XlO 10 

09 0+2.4 

04.o_ 79 

(3.46+??) xlO 9 

- 

- 

LBG J095956.93+022118.5 

2.470 

UV 

(2.48/g'g 3 ) X 10 10 

IS.?/®;® 

(1.77+??) xlO 9 

- 

- 

LBG J100059.45+021957.4 

2.470 

Hcn/Lya 

(3.75/g'gg)Xl0 10 

16-9 

(2.66+??) XlO 10 

- 

Opt/X-ray 

LBG J100100.91+021728.1 

2.470 

Lya 

(7.18++)xl0 10 

30.71+ 

(2.84+1!) xlO 9 

- 

- 

DSFG J100018.17+022250.4 

2.470 

Ha 

(l-Ol/ggsJXfO 11 

365+?? 

(1.91++)xl0 10 

Merg 

- 

LBG J100014.24+022516.7 

2.471 

UV 

(3.24++) xlO 10 

46.711;® 

(1.13+'??)xl0 10 

Sph/lnt 

- 

LBG J100115.18+022349.7 

2.471 

Lya/UV 

(3.48lg5g)Xl0 10 

17.8+;? 

(5.54++)xl0 10 

- 

Radio 

DSFG J100016.57+022638.4 

2.472 

Ha 

(8.82++) xlO 9 

714+ 210 

(1.09++)xl0 n 

Disk 

Radio 

DSFG J100056.83+022013.3 

2.472 

Ha 

(3.7l/g'gg)Xl0 10 

1961®? 

(8.13+'?|)xl0 10 

- 

- 

LBG J100018.04+021808.6 

2.472 

UV 

(4.58+??) XlO 10 

71 7+8-5 
11,1 - 17.1 

(7.44+-??) xlO 9 

Sph/lnt 

- 

LBG J100020.50+022421.5 

2.472 

UV 

(4.37++) XlO 10 

33 . 1 I + 3 

(1.94+1?) XlO 10 

Sph 

- 

LBG J100109.29+022221.5 

2.473 

Lya/UV 

(3.73/g'gg)Xl0 10 

52.6l+i 

(1.04+??) xlO 10 

- 

- 

LBG J100015.38+022448.3 

2.474 

Lya/[Civ] 

(4.93+?!) XlO 10 

11 . 81 ?;? 

(1.32+'??)xl0 10 

Sph 

UV 

LBG J100033.20+022225.0 

2.474 

Lya/UV 

(2.48++) XlO 10 

63.2li 5 9 6 

(1.52+?!)xl0 10 

Sph 

- 

LBG J100008.88+023044.1 

2.475 

UV 

(4.45/g'gg) X 10 10 

16.2+? 

(1.93++)xl0 10 

- 

- 

LBG J100012.37+023707.6 

2.475 

UV 

(3-05++) XlO 10 

I 31 I 32 

(2.16++)xl0 10 

- 

- 

LBG J100015.87+021939.5 

2.475 

UV 

(3.31++)xl0 10 

79-411+ 

(3.15+'??)xl0 10 

Disk/Int 

- 

LBG J100025.28+022643.3 

2.475 

UV 

(4.58/g'gg)Xl0 10 

45.2+;? 

( 2 . 01 ++)xl 0 10 

Merg 

- 

LBG J100116.15+021854.2 

2.475 

UV 

(3.86/g'gg)Xl0 10 

90 o + l -3 
— 12.5 

(4.99+??) xlO 9 

- 

- 

LBG J100002.03+023012.9 

2.477 

Lya/UV 

(5.46+??) XlO 10 

90 q+ 1-9 

*o.o _ 14 Q 

(4.38++) xlO 9 

- 

- 

LBG J100054.07+022104.4 

2.478 

Lya/[Civ] 

(4.37+??) XlO 10 

22.91+ 

(2.00+'!?) xlO 10 

- 

UV/X-ray 

DSFG J100026.73+022411.3 

2.478 

Ha 

(7.55++)xl0 9 

284+™ 

(1.64++)xl0 10 

Disk 

Opt 

LBG J100024.21+022741.3 

2.479 

UV 

(5.47+'??) XlO 10 

i 8 .ol?;? 

(1.04+!!) xlO 10 

Sph 

- 

DSFG J100018.21+023456.7 

2.479 

CO(l-O) 

<3.Ox 10 ll) 

4001J+ 

(2.13+'??)xlO n 

- 

- 

QSO J100021.96+022356.7 

2.480 

Lya 

(2.62++X10 11 


(6.59+'??)xl0 10 

QSO/Sph UV/X-ray 

DSFG J100027.14+023140.8 

2.483 

Ha 

(6.04/g/g)Xl0 10 

1090+?? 

(4.10+'??)xl0 10 

Sph 

Opt 

LBG J100004.33+022654.1 

2.480 

UV 

(6.22+'??) x 10 10 

37 9“*" 3 - 6 

01 ‘ v - 8.3 

(4.03+??) XlO 9 

- 

- 

LBG J100033.91+022713.2 

2.481 

UV 

(2.17/g'gg)Xl0 10 

27.3+;? 

(2.57+1?) xlO 9 

- 

- 

LBG J100023.31+023537.5 

2.487 

UV 

(3.69/°®®) xio 10 

27 9+ 0,3 

(3.96+?!) XlO 9 

- 

- 


Source names indicate the selection method: submillimeter selection by DSFG, optical by LBG, and one unequivocal quasar 
by QSO. Luv is the rest-frame ultraviolet luminosity computed from an extrapolated apparent magnitude at rest-frame 1600A. 
Morph Class refers to the CANDELS visual classification scheme applied to HST H-band imaging. ’Sph’ refers to spheroids, 
’Disk’ refers to disk-like galaxies, ‘QSO’ refers to an unresolved point source, ‘Merg’ refers to a galaxy merger, and ’Int’ refers 
to some signature of galaxy interactions taking place. 


out, galaxy cutouts for cluster members and control were 
scrambled to ensure unbiased results. Although limited 
by small numbers, we find that 10 of 21 proto-cluster 
member galaxies (48±10%) appear to be irregular or un¬ 
dergoing interaction while only 5 of 25 control galaxies 
exhibit interaction (20±8%). Even with removal of the 
DSFGs, a high interaction fraction (7/16=44%), is found 
for the proto-cluster members. 

3.5. Estimating Halo Mass. 

To estimate the net dark matter halo mass of this 
structure, we use abundance matching techniques from 


large-volume simulations (Behroozi, Wechsler, & Con-| 
roy 2013). Due to PCL1002’s filamentary structure, 


filamentary 


we estimate the total halo mass by summing the es¬ 
timated halo masses for each galaxy in the structure 
using its stellar mass. We estimate a lower limit of 
Mhaio >(8±3)xl0 13 M 0 at z = 2.47. Using a model 
for mass-dependent exponential growth, we estimate the 
z = 0 halo mass to be (2± l)xl0 15 M 0 ( Wechsler et al.| 


2002 Chiang et al. 2013), about twice as large as the 
Coma supercluster (which has a mass of ~lxl0 15 Mg). 
Note that this dark matter halo mass estimation method 
differs from others which assume linear bias and splier- 
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ical collapse models ( Mo fe White|[l996 Peacock 19991 


which we suggest are not applicable to z Xl 2 filamentary 
structures. 


Miller et al. 

2015) 

s in DSFGs c 

o not 


trace massive proto-clusters. They argue that such sig¬ 
nificant over-densities are due to “merger bias” whereby 
the submillimeter-luminous phase is too short-lived and 
rare to be a useful probe of the most massive halos at 
high- 2 . Indeed, there are several massive proto-clusters 
at 2 > 2 that contain no DSFGs, and some structures 


DSFGs (iCapak et al. 

2011’; Hodg 

e et al. 

2012 

Wal- 

ter et al. 

2012)'.' The 

Miller et al. 

(20151 work argues 


Fig. 3.— 4" X 4” rest-frame optical cutouts for PCL1002 members from HST- WFC3. Galaxy names are indicated along with morphological 
and interaction indicators (see Table |2|. 

rare. Even scaling density estimates proportional to the 
LBG overenhancement, the DSFG/AGN presence is ex¬ 
cessive. The extra DSFG/AGN can be explained with 
only two possible physical interpretations: either the DS¬ 
FGs and AGN are short-lived and triggered simultane¬ 
ously via a process related to the over-dense environment 
or the DSFGs and AGN must be much longer lived than 
expected by existing observations of similar phenomena 
in the field. 

The evidence we present here is suggestive of the for¬ 
mer, that DSFGs and AGN are short-lived, and in that 
case, their triggering must be environmentally driven. 
This is supported by the increased interaction frac¬ 
tion seen in proto-cluster members’ optical morpholo¬ 
gies, larger ISM masses in proto-cluster members (al¬ 
beit a marginal distinction requiring confirmation), and 
the lack of physical motivation for long Q SO lifetimes 
(Martini 2004 Hopkins & Hernquist [2009 ) beyond the 
increased gas-feeding argument often argued for DSFGs. 
If correct, these observations provide the first concrete 
evidence that environmental triggering can occur over 
extraordinarily large volumes (15000 Mpc 3 comoving) at 
2 > 2. 

It should be noted that the structure’s member LBGs 
do not appear to have sufficiently different observational 
characteristics than LBGs in the field. Comparison 
against the control sample indicates that proto-cluster 
LBGs have slightly higher stellar masses (by a factor 
of 1.5±0.3) and similar star-formation rates. Neverthe¬ 
less, as a sheer consequence of their number, LBGs dom¬ 
inate the calculation of the structure’s net dark mat¬ 
ter halo mass which is estimated to be significant at 
>8xl0 13 Mg. This is predicted to mature to a few 
xlO 15 M 0 at 2 ~ 0. 

Cosmological dark-matter simulations suggest that to¬ 
day’s most massive c lusters occupied vol umes several 
hundred times larger (Onorbe et al. 2014) at 2 ~ 2.5, 
having not yet virialized into the compact structures 
we see today. PCL1002 affirms these predictions, as its 
15000 Mpc 3 volume is predicted to collapse to a volume 
of « 50 Mpc 3 at 2 ~ 0, which is consistent with all nearby 


massive over-densities at high-redshift because Poisson 
noise dominates for low numbers of DSFGs. While this 
is likely the case for 1-3 DSFGs, our results (and our 
interpretation of their Figure 3) imply that the oppo¬ 
site is actually true for sufficiently large numbers of DS¬ 
FGs per proto-cluster. Instead we suggest that spec¬ 
troscopic incompleteness in both LBG and submillime¬ 
ter samples around high -2 proto-clusters has potentially 
hindered the discovery of more starburst-enriched and 
AGN-enriched proto-clusters. If this is the case, aggres¬ 
sive spectroscopic follow-up of DSFGs might substan¬ 
tially help the effort to identify high -2 proto-cluster en¬ 
vironments, where the spectroscopic investment needed 
for LBG samples is prohibitive on large scales (i.e. iden¬ 
tification of >40 spectroscopically-identified LBGs, with 
a resulting <5lbg ~ 3 — 5). With the identification of four 
20-30 arcminute-scale over-densities containing >6 DS¬ 
FGs (HDF, PCL1002 and SSA22) in only a few degrees 
of deep submillimeter surveys, the potential to identify 
further massive cluster progenitors via their member DS¬ 
FGs and AGN is quite promising. 

4. DISCUSSION & CONCLUSIONS 

The over-abundance of what are thought to be short- 
duration, --1 00 Myi'. phenomena (both DSFGs and lu- 
minous AGN; Bothwell et al.||2013[ Marconi et al.|[2004 1 
in an extended proto-cluster structure is conspicuously 
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> 5 x 10 14 Mq clusters. While this confirms the notion 
of a genuine massive cluster in formation, this volume¬ 
scaling also implies that most clusters will have similarly 
large sizes at 2 > 2.5, subtending areas half a degree 
across on the sky. Observationally isolating massive clus¬ 
ters in formation then comes down to accurate redshift 
identification to within A z = 0.03, because other obser¬ 
vational characteristics are not sufficiently environmen¬ 
tally distinct at these epochs. 

Identifying and correctly classifying PCL1002 as a 
massive galaxy cluster progenitor would not have been 
possible without the concentrated presence of cosmologi- 
cally rare phenomena like dusty starbursts and luminous 
AGN. Future large and deep £ 100 deg 2 submillimeter 
surveys could play a crucial part in statistically charac¬ 
terizing the population of such large structures during 
their formation epoch, as groups of DSFGs and lumi¬ 
nous AGN can act as signposts to the largest mass con¬ 
centrations. Equally important is complete spectroscopic 
follow-up at z > 2 over those wide-a reas, like large opti¬ 
cal campaigns similar to HETDEX (Chiang et al. [2014 1 
and potential future large millimeter line searches tar¬ 
geting CO or [Cil] with a wide-bandwidth spectrometer. 
The discovery of more high- z , starbursting over-densities 
will be inevitable, but it will be the constraints on their 
volume density and assembly timescale that will have 


significant repercussions on cosmological hydrodynamic 
simulations and the formation mechanisms of the Coma¬ 
like superclusters we see today. 
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